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I. INTRODUCTION 

Since the experimental synthesis of graphene the inter- 
est in this material has mostly focussed on its electronic 
properties As time goes by and as the experiments 
become more accurate the interest is being displaced to 
the morphological aspects and their possible influence on 
the electronic properties . Part of the interest on these 
issues is prompted by the necessity to find the mecha- 
nisms limiting the mobility of the graphene samples 3], 
an important issue for the potential technological appli- 
cations. Another source of interest lies on the special 
mechanical and clastic properties of graphene recently 
tested in experiments [4, 5, 6, 7, 8] and on the fascinating 
possibility to explore the physics of real two dimensional 
crystals. 

Local disorder in graphene has been thoroughly stud- 
ied and we refer to the review article [l[ for a fairly com- 
plete list of references. A different type of disorder is pro- 
vided by the observation of ripples in suspended grap hene 
and in graphene grown on a substrate [Til [l^ . 
The elastic and mechanical properties of graphitic struc- 
tures have been studied intensely in the past, mostly in 
the context of understanding the formation of fuUerenes 
and nanotubes. Very little work has been done for the 
flat graphene sheet [H, and topological defects have 
been often excluded in these studies. In the fuUerene 
literature it was established that the formation of topo- 
logical defects (substitution of a hexagonal ring by other 
polygons) is the natural way in which the graphitic net 
heals vacancies and other damages produced for instance 
by irradiation [l^. Among those, disclinations (isolated 
pentagon or heptagon rings), dislocations (pentagon- 
heptagon pairs) and Stone- Wales (SW) defects (special 
dislocation dipoles) were found to have the least forma- 
tion energy and activation barriers. Dislocations and 
SW defects have been observed in carbon structures 
and are known to have a strong influence on the elec- 
tronic properties of nanotubes. Dislocations in irradiated 
graphitic structures have been described in [l7[ and ex- 
perimental observations have been reported in graphene 



grown on Ir in fisll . 

In previous works we used a general relativity formal- 
ism to study the influence of curved portions of the lattice 
on the electronic properties of graphene by coupling the 
Dirac fermions to the corresponding curved space. We 
found that curvature generates a fictitious gauge field 
- a result that is also obtained in more conventional 
descriptions - and, in addition, we predicted a space- 
dependent parameter in the kinetic term (similar to a 
space-dependent Fermi velocity) that depends on the in- 
trinsic curvature of the sarnple and has non-trivial effects 
on the density of states [H, [13, Hi], [l^l and on the conduc- 
tivity j23| . In the present work we extend this formalism 
to include a finite density of dislocations. Dislocations 
are modelled by adding torsion to the graphene sheet 
either curved or flat. 

Torsion was first considered by Cartan and then by 
Sciama and Kibble in order to deal with spin in General 
Relativity [l^l . Within this scheme the spin of a particle 
turns out to be related to the torsion just as its mass is 
related to the curvature. Torsion is generally ignored in 
the context of GR, because the torsion field is algebraic in 
the Einstein-Cartan theory (its equations of motion are 
constraints, it has no dynamics) and as a consequence 
there is no torsion in the absence of matter. When mat- 
ter generating torsion is added, its observable effects are 
suppressed by the smallness of the gravitational coupling 
and they are negligible These two drawbacks are 

overcome when torsion is included in a geometric fashion 
as associated to dislocations. In the context of electron 
systems like graphene the dynamics of the lattice defects 
occurs at a much higher energy than the electronic pro- 
cesses (in graphene, lattice processes are related to the 
sigma bonds and have typical energies of the order of tens 
of eV while the continuum relevant low energy processes 
are of few tens of mcV). It then makes perfect sense to 
consider the motion of electrons in a frozen geometry 
and hence the fact that torsion has no dynamics is not 
a problem. Moreover the suppression of the coupling of 
the torsion and electronics degrees of freedom will be of 
the order of magnitude of the quotient between electronic 
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and elastic degrees of freedom that, although small, will 
be much bigger than the gravitational scale. 

The purpose of these notes is twofold: First we describe 
a continuum description that allows to explore the effects 
of dislocations on the electronic properties of graphene 
and other materials, and, second, we propose an "analog" 
model that amplifies the effect of torsion in a General 
Relativity context. 

II. DISLOCATIONS AND TORSION 

Topological defects in crystal lattices can be classified 
in two kinds: Dislocations and disclinations. A disloca- 
tion can be thought of being formed by performing a cut 
in a bulk material and introducing extra rows of atoms 
(Volterra process). These types of point defects can be 
easily identified by going around a closed circuit enclos- 
ing the end of the line cut with discrete lattice steps. 
Without defects the path will end up at the same point. 
In the presence of a dislocation it ends up somewhere 
else and an extra vector has to be added to close the 
path. This is called the Burgers vector. Dislocations can 
be further classified by the direction of the Burgers vec- 
tor with respect to the plane defined by the circuit: In 
an edge dislocation the Burgers vector lies in this plane, 
while in screw dislocation it is perpendicular to it. Fig. 
[1] shows an edge a and screw dislocation in the square 
lattice characterized by a Burgers vector parallel (per- 
pendicular) to the displacement. A disclination is simi- 
larly formed by adding a wedge of atoms, and it can be 
identified by computing the total angle subtended by a 
closed path enclosing the defect (Frank angle). A discli- 
nation dipole can also be seen as a dislocation. ^From the 
point of view of geometry, dislocations represent a source 
of translational mismatch "anholonomy" , (see Appendix 
I A 21 for a rigorous description of torsion), and as we will 
see, are very naturally described in the continuum by a 
connection with torsion. In differential geometry, a con- 
nection defines the notion of parallel transport. When a 
vector BP is parallel transported from to + dx^, it 
experiences an infinitesimal variation given by 

SBP = -TPp{x)B''dxP. (1) 

The equivalent of the closed circuit in the definition of the 
Burgers vector can be defined by taking two infinitesimal 
vectors mP and n'^ , and trying to build a parallelogram 
with them. For this purpose, one parallel transports the 
vector niP along n'^: 

m'P + nP ^TP^^mPn", (2) 

and the vector n'^ along mP 

n'P ^mP + nP -TP^mPn". (3) 

As depicted in Fig. [Ill if there is a dislocation in the 
region enclosed by the path, the parallelogram obtained 




FIG. 1: (Color online) Upper: Edge dislocation in a square 
two-dimensional lattice. The Burgers vector is parallel to the 
displacement. Down: Screw dislocation in a three dimen- 
sional cubic lattice. The Burgers vector is perpendicular to 
the displacement. 



does not close, and the part that is missing is propor- 
tional to the antisymmetric part of the connection, de- 
fined as the torsion: 

m'P - n'P = [r^, - re J m^n'^ ee T^^m^n'^ (4) 

Notice that the two paths would have closed in a curved 
space with a standard affine connection. 

In the continuum limit, one can describe a density of 
dislocations by the density of Burgers vector, and we see 
that torsion plays precisely this role. Hence the natural 
identification of torsion with dislocations. The connec- 




FIG. 2: Non-closure of infinitesimal parallelograms due to the 
presence of torsion 



tion of torsion with the continuum theory of crystal dislo- 
cations goes back to Kondo [26| and has been formalized 
in [23, [l!]. Various aspects of the problem have been 
explored in [2^ . A nice review on the relation of gravity 
with topological defects in solids is [ 3(]|| . The geometric 
approach to defects in solids [1^, l3lj| relates the metric 
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FIG. 3: (Color online) Structure of the glide (up) and shuffle 
(down) dislocations in the planar graphene lattice. 



of the curved -crystalline- surface with the deformation 
tensor and establishes that disclinations are associated to 
the curvature tensor and a finite density of dislocations 
generates a torsion term. 

In flat surfaces screw dislocations do not exist. Dislo- 
cations in graphene are made of pentagon-heptagon pairs 
and they have been widely studied in connection with the 
properties of carbon nanotubes ^] and, more recently, 
in the flat graphene sheets or ribbons [H, [s^, • Ob- 
servations of topological defects in graphene have been 
reported in |16, 36, 37]. The stability of dislocations in 
graphene has been explored recently in [s^l and two types 
of stable dislocations were found. The most common 
are the so-called glide dislocations made by a pentagon- 
heptagon pair shown in the upper part of Fig. [Ill In ad- 
dition we can have shuffle dislocations (lower part of Fig. 
im made by a octagon with a danghng bond. Both are 
edge dislocations and add one or several lines of atoms to 
the lattice. The difference is that shuffle dislocations con- 
tain an extra atom with respect to the glide dislocations. 
In an infrared continuum formulation both are seen as 
point defects that do not disturb the lattice out of the 
core. The difference with conical defects as pentagons or 
heptagons is that when an electron circles such a defect 
it will found itself displaced from the original position a 
distance proportional to the Burgers vector as shown in 
Fig. M 



III. COUPLING DIRAC FERMIONS TO A 
CURVED SPACE WITH DISLOCATIONS 

The model that we propose to study the low energy 
properties of graphite and graphene in the presence of 
lattice defects (disclinations and dislocations) is based 
on the assumption that the low-energy description of 



the quasiparticles in terms of massless Dirac fermions 
remains valid in the presence of a diluted distribution of 
defects. The stability of the Fermi points in multilayer 
arrangements of the honeycomb lattice has been demon- 
strated in [s^l. As we have described in the previous 
section the presence of dislocations in a crystal can be 
described in the continuum limit by adding torsion to the 
space spanned by the lattice. We will here follow a covari- 
ant approach close to the general relativity formulation 
what allows to consistently couple matter to geometry. 
This description began with the early paper of ref. [39| 
where the Dirac equation on a sphere was used to solve 
the electronic spectrum of fuUerenes. The pentagonal de- 
fects (disclinations) were described by a fictitious gauge 
field. The formalism was latergeneralized to study the 
influence of disclinations [iO, |4l[ and of smooth curva- 
ture [111, [l^l on the density of states of graphene. These 
previous works showed that the covariant formulation al- 
lows to obtain sensible predictions that are in agreement 
with the physical observations when the infrared proper- 
ties are considered. In particular it was shown that the 
curvature induces a fictitious vector gauge field. We will 
here analyze the nature of the effective gauge fields cre- 
ated by the presence of dislocations within the covariant 
framework. 

Before introducing the changes induced by having tor- 
sion wc summarize the situation in a curved space. 

The dynamics of a massless Dirac spinor in a curved 
spacetime is governed by the modified Dirac equation 
(see Appendix E]): 



z7^(r)^^V = 



(5) 



where 7'^ are the curved space 7 matrices that depend 
on the point of the space and can be computed from the 
generalized anticommutation relations 

{7'^(r),7'^(r)} = 25'^''(r), 

and the covariant derivative operator is defined as 



D„ 



da - 



where £7^ is the spin connection of the spinor field whose 
physical implications have been discussed in [21I, [13] ■ 

Coupling a Dirac field to a curved space with torsion 
gives the so-called Einstein-Cartan-Dirac theory. The de- 
tails of the derivation of the Dirac equation in a space 
with torsion are given in Appendix lA 41 Here we will 
simplify the description putting it in a different, more 
physical way. The minimal coupling of any geometrical 
or physical fields to the Dirac spinors adopts always the 
form of a covariant derivative: 



da 



Da 



da + Aa 



(6) 



where the given vector can be an electromagnetic poten- 
tial induced by a real electromagnetic field or any other 
real or fictitious gauge field associated to deformations 
or to geometrical factors 42|. In the case of having a 
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density of dislocations in the graphene sheet modelled 
by torsion we can construct two potential candidates to 
gauge fields. In four dimensions, from the rank three 
torsion tensor (jA14[) T^,yp the following vectors can be 
built: 

^fi 9 ^^upp 1 S ^ pu paTf per ■ ('^) 

The field is a real vector and can be associated to the 
density of edge dislocations while Sp is an axial (pseudo) 
vector associated to the density of screw dislocations. 
These fields couple to the vector and axial current density 
respectively so the full Lagrangian is: 

iint = j d^x^ [^{d^ + leV^ + iiii^'S^)]^, (8) 

where e and rj are coupling constants related to the den- 
sity of edge and screw dislocations respectively. The 
physical consequences of the vector torsion are then sim- 
ilar to these coming from curvature or elasticity. The 
coupling of the spinor to the axial field can have im- 
portant consequences as it breaks time reversal symme- 
try. 

A. Two dimensional systems 

In (2-1-1) dimensions the completely antisymmetric 
part of the torsion tensor is proportional to the Levi- 
Civita tensor: 

^pvp — ^pvp^: (9) 

where $ is a pseudoscalar field. 

In 2-1-1 dimensions there are no screw dislocations, sim- 
ply because there is no third spatial dimension, and the 
antisymmetric part of the torsion is rather related to time 
displacements. Since the time components of the torsion 
generated by a dislocation are zero, this term vanishes 
always. The remaining coupling to the trace part acts a 
new vector gauge field with special symmetry properties 
associated to the dislocations structure of the sample. 
This field will add up to these coming from curvature or 
to real electromagnetic fields. 

A remaining interesting question is that of how this 
extra gauge field couples to the electronic excitations 
around different Fermi points (valleys). It is known that 
the fictitious gauge fields arising from elastic deforma- 
tions as well as these coming from curvature, give rise to 
effective magnetic fields pointing in opposite directions 
in different valleys. The reason is that the sign of the 
coupling of the effective vector fields depends uniquely 
on the definition of the spin connection (jA6p 

np^lio;''hanb], (10) 

and is determined by the product of gamma matrices 
in it. As the effective Hamiltonians around each of the 



Fermi points differ in the sign of only one of the gamma 
matrices (they are similar to parity conjugated) the spin 
connection has opposite signs in the two valleys. This 
is irrespective of whether the connection has or not an 
antisymmetric part and hence applies also to the coupling 
induced by torsion. 



B. Three dimensional systems supporting Dirac 
fermions 

We have seen in the previous sections that the effects of 
dislocations in two dimensional crystals like graphene are 
- at most - describable in terms of a fictitious gauge field 
similar to the one associated to curvature. The formalism 
outlined in this work acquires full power when applied to 
three dimensional crystals supporting screw dislocations 
and having Dirac fermions in the spectrum. Graphite 
is one of the most obvious candidates where screw dis- 
locations are very common [isj and elementary excita- 
tions behaving as Dirac fermions were predicted near the 
K points of the band structure in the early calculations 
p4 , |45| . These Dirac fermions have been experimentally 
confirmed in [i^, \^ . 

The other perfect candidates are the the topological 
insulators [48l.l49l| specially the Z2 type that are protected 
by time reversal symmetry. The axial part of the torsion 
will couple to the three dimensional Dirac fermions, break 
the TRS and affect the topological stability. The density 
of states associated to the screw dislocations in [s^ can 
be reproduced with the techniques outlined in this work 
and will be reported elsewhere [6l| . 



IV. DISLOCATIONS IN CONVENTIONAL 
ELECTRONIC SYSTEMS 

Another interesting question is how dislocations affect 
the electronic properties of crystal systems whose quasi- 
particles are described by the conventional Schrodinger 
equation. Within the present formalism the question can 
be posed as: Do scalar fields see torsion? 

The question is interesting and relevant since the ma- 
jority of electronic systems (Fermi liquids) obey a con- 
ventional Schrodinger equation and can be modelled with 
scalar fields. The point to notice here is that, although 
scalar fields transform in quite a trivial manner under ro- 
tations and translations, building a proper Hamiltonian 
means building derivatives and there is where curvature 
and perhaps torsion plays a role. The Lagrangian for a 
scalar field in a curved space is 

Co = J (11) 

As we see it involves only the metric and hence it does not 
see the torsion. In the spirit of building phcnomcnological 
couplings of the electronic density to the dislocations the 
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most natural scalar couplings that we can form from the 
torsion vectors ([7]) are 



of the effective gauge fields induced by dislocations on 
graphene have been discussed in (H^, [s^, [H, [13] ■ 



(12) 



where is the vector associated to the trace of the tor- 
sion, i. e., to the density of edge dislocations. The cou- 
pling §2 being quadratic in the torsion will be suppressed 
with respect to gi by an extra power of the density of 
dislocations. In the two dimensional case we can form 
the same interactions with the scalar field of 

Landau levels in the presence of a screw dislocation 
have been explored with a similar formalism in [2^ . 



V. CONCLUSIONS AND DISCUSSION 

The proposal of modelling dislocations in a crystal with 
a geometrical approach including torsion within the elas- 
ticity theory is very old. What is new in our approach is 
the suggestion to study the influence of dislocations on 
the electronic properties of graphene and other electronic 
systems having Dirac fermions in the spectrum, by cou- 
pling the Dirac equation to a generalized curved space 
with torsion. In what two dimensional graphene sheets 
are concerned the main results of this work are somehow 
disappointing: spinors in two spatial dimensions do not 
couple to the axial part of the torsion (since there are 
no screw dislocations). Such a coupling would have had 
very important consequences. The coupling to the trace 
part is in the form of a vector gauge field similar to the 
one produced by the curvature or by strain in a more 
conventional approach. The spacial distribution of this 
new gauge field will nevertheless depend on the density 
of dislocations. 

The most interesting applications of the formalism out- 
lined here will be in the area of topological insulators that 
often support massless chiral fermions in three dimen- 
sions. In ^51] it was shown that the coupling of spinors 
to the axial vector induced by torsion generates at higher 
order a four Fermi interaction on the fermions what had 
an influence on neutrino oscillations. In the context of 
condensed matter such a coupling would renormalize the 
Hubbard interaction or generate it if it was absent. The 
presence of the torsion axial vector in the Dirac equa- 
tion will have also implications on various fundamental 
aspects of the physics of condensed matter system de- 
scribed by Dirac fermions. As it breaks time reversal 
symmetry, it can prevent weak localization. As for the 
conductivity, the disorder induced by a random distribu- 
tion of dislocations may change the universality class of 
previous models based on vacancies or other T-preserved 
disorder. The breakdown of time reversal symmetry will 
also lift the quantum protectorate associated to the Fermi 
points [111 and allow the opening of a gap although the 
torsion itself will not explain the observations reported 
on gaps in graphene (s^ . [53| . The physical implications 



APPENDIX A: THE DIRAC EQUATION IN A 
SPACE WITH TORSION AND CURVATURE 

1. Dirac equation in a curved space without torsion 

The behavior of spinors in curved spaces is more com- 
plicated than that of scalar or vector fields because their 
Lorentz transformation rules do not generalize easily to 
arbitrary coordinate systems Instead of the usual 

metric g^j/ we must introduce at each point X described 
in arbitrary coordinates, a set of locally inertial coordi- 
nates the vielbein fields e°(a;), a set of orthonor- 
mal vectors labelled by a that fixes the transformation 
between the local and the general coordinates: 



=x- 



(Al) 



The curved space gamma matrices 7^(x) satisfying the 
commutation relations 

{7^7,} = 2g^,, (A2) 

are related with the constant, flat space matrices 7° by 

7^(a;) - e^^7^ (A3) 

The spin connection 17^ (x) is defined from the fielbein by 

1 



(A4) 



where is the covariant derivative acting on the 



vectors as 



pA a 



(A5) 



The spin connection for 1/2 spinors can be written as 
1 



-UJ 



ah 



[7a, 76] 



where u^'' are the spin connection coefficients 



ah 



(A6) 



(A7) 



In the flat case r„„ is the usual affine connection 



(Christoffel symbols) related to the metric tensor by 



1 



dgo 



dgtiu dg^ 



fjLCT 



}• 



(A8) 



2" ' dx^" dx" dx" 
Finally, the determinant of the metric needed to define a 
scalar density lagrangian is given by 

^/=^=[det(g^.)]l/2 ^ Aet[el{x)]. (A9) 

These formulas apply to a curved space without torsion. 
The Dirac equation is then written as 

i7^i?^V = 0, (AlO) 
where = 9^ + is the covariant derivative acting 



on the spinors in the curved space and 7^ 
curved gamma matrices. 



7''e^ the 
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2. Non-singular torsion in differential geometry 

In differential geometry we consider a curved manifold, 
described by a set of coordinates x^^. To this manifold 
we associate a metric g^j^^, which defines the distance lo- 
cally. A notion of parallel transport is also needed to 
compare vectors in different points, in order to construct 
meaningful derivatives. The connection and the metric 
are independent geometric objects that are usually re- 
lated by the so-called metricity condition ensuring that 
the relative angles and magnitudes of vectors are pre- 
served under parallel transport. The metricity condition 
reads: 



r^.ffPM = 0, (All) 



this is, the covariant derivative of the metric vanishes. 
The metricity condition fixes the symmetric part of the 
connection and leaves the antisymmetric part undeter- 
mined. The most general connection, compatible with 
the Riemannian metric, is 



A 



K 



(A12) 



A 



where: j- are the Christoffel symbols and K^^^ is 

called the contortion tensor related with the torsion ten- 
sor by 



The torsion is the antisymmetric part of the connec- 
tion: 



(A14) 



This means that when choosing a connection for a met- 
ric space, we are always free to choose the antisymmet- 
ric part of the connection (the torsion) at will. Only 
the symmetric part is fixed by the metricity condition. 
The usual choice is to take the torsion equal to zero, so 
that the metric alone determines all geometric proper- 
ties. This is, for example, the case of General Relativity. 
The only restriction on K is that 

Kx^u = -K,px (A15) 

hence in a d-dimensional space it has ^[d?'{d — 1)] com- 
ponents. The differential geometry of a general curved 
space is determined by the two tensor fields g^^, and T'j^^ 
(or, equivalently if^j,). 

In a space with torsion the shortest path (geodesies) 
and the covariantly paralell path do not coincide and, as 
was mentioned in Section |TT1 parallelograms do not close. 



3. Torsion in a 2D flat space 

Aiming to a description of dislocations alone, we study 
here the special case of a space with torsion but no cur- 
vature, in 2-1-0 dimensions. This will be the spatial part 



of the space-time used to describe graphene with dis- 
locations as the Dirac equation in a space with spatial 
torsion. 

As we will see, life is very simple in this space. First 
of all, since there is no curvature, there is a global coor- 
dinate system in which the metric tensor is the identity. 
Or seen from the tetrad perspective, there is a coordi- 
nate system in which the tetrads are trivial, and thus 
the Lorentz connection and the Einstein connection are 
equal. (Since the tetrads are trivial, there is no difference 
between latin and greek indices now, and moreover, since 
we are in 2-|-0, upper and lower indices are equal.) In this 
particular case, with the tetrads chosen in this way, we 
can say that the torsion in encoded in the connection. 

Moreover, the antisymmetry of the torsion tensor in 
2D immediately implies that it can be written as: 



(A16) 



where is an arbitrary vector field, basically the density 
of Burgers vector. It could be thought that this vector 
has some further constraints due to the Bianchi identities 

identities are trivial 



but it can be checked that in 2D if Rp,^pa = these 



4. Dirac equation in a curved space with torsion 



As have seen in appendix lA 21 inclusion of torsion in 
the geometrical description simply changes the connec- 
tion and hence the covariant derivative. 

It would be tempting then to propose that the Dirac 
equation in the space with torsion and curvature is given 
by eq. ([5]) with the appropriate connection but there are 
some important subtleties that we will specify [1^ . 

Consider the Dirac field in flat space. The massless, 
manifestly hermitian Dirac Langrangian may be written 
as: 



Noting that: 



(A17) 



(A18) 



plus the fact that a total derivative doesn't affect the 
equations of motion (with suitable boundary conditions), 
the Lagrangian (|A17[) can be traded for: 



C 



I 



(A19) 



Indeed, the Dirac equation can be derived by variation 
with respect to tp of any of those Lagrangians. The sim- 
pler form (|A19P seems to be not hermitian, but this is not 
a problem since we can always introduce the boundary 
term back. 
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Now consider the same problem in a manifold with cur- 
vature and torsion. Consider the hcrmitian Lagrangian: 

C^J d"x^^ [^^D^i; - (D^^p) , (A20) 

with = ~ ^fj.bcl^l'^- Now the equivalent of the 
chain rule is: 

(D^^P) 7^7A = [^^^^] - ,PD^ [Y^] . (A21) 

Note the following important point: since the covari- 
ant derivative contains gamma matrices, we can not just 
commute it to make it act directly on the spinor field. 
The "chain rule" that we have to use is really: 

(A22) 

where we have used: 

^b^c^a ^ ^a^b^c ^ 2 (^\<^<^ - ^'^.fb^ . (A23) 

The commutation with the gamma matrices has in- 
troduced a new term, the trace of the connection. This 
means that the equivalence of the simpler form given in 
eq. (jA19p is not just obtained by promoting the deriva- 
tive to a covariant one, but it also requires the introduc- 
tion of this trace [5§| . 

In four dimensions, the following identity: 

^a^b^c ^ ^agbc ^ ^Cgab _ ^b ^ac _^ i^abcd^5^d^ (^24) 

and the definition of the connection (|A6p allows to rewrite 
this as: 

£ = J d^x (^i^j'^daiJ + ^V^t'^tS'V^ , (A25) 

or: 

£ = fd^x {^rdai^ + te^'^^n^di^^^'ri^) , (A26) 



Which reveals the well known result that, in General Rel- 
ativity, fermions only couple to the antisymmetric part of 
the torsion [sl]. In three dimensions using the identity: 

[7",7''] = -2e'^Nc, (A27) 
we can write the action in the form 

£ = y" d^x (i^7"9aV' + ^e^'^Tahc^) , (A28) 

from where the Dirac equation can be extracted directly 
to read 

[^da + ze°'^Tafc,]V^ = 0, (A29) 



Note added: After completing this work we became 
aware of ref. ^6d\ where similar issues are discussed. We 
think that, although there is some obvious overlapping, 
the two works entertain different points of view and can 
be seen as complementary. 
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